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Abstract. We address the point-to-face approximate shortest path prob- 
lem in Ti?: Given a set of polyhedral obstacles with a total of n vertices, 
a source point s, an obstacle face /, and a real positive parameter e, com- 
pute a path from s to f that avoids the interior of the obstacles and has 
length at most (1 + e) times the length of the shortest obstacle avoiding 
path from s to /. We present three approximation algorithms that take 
0(n*(L-Hlog(n/e))Ve2 _^n2(L + log(n/e))/e^) time, 0(rp_p(n) * (l/e^)) 
time, and 0(n^A(n) log(n/e)/e'' -I- log(np) log(n log p)) time, respec- 
tively, where L is the precision of the integers used, Tp-p{n) is the time 
complexity of the point-to-point shortest path algorithm used, p is the 
ratio of the length of the longest obstacle edge to the Euclidean distance 
between s and /, and A(n) is a very slowly-growing function related to 
the inverse of the Ackermann's function. 



1 Introduction 

The Euclidean shortest path problem among obstacles in the plane or space is 
one of the oldest and well-known problems in computational geometry. It has 
been intensively studied, see pilS HnifTSHroil^l^ . as well as the survey by 
Mitchell [11]. In general the problem is stated as: Given a set of obstacles in the 
d-dimensional space TZ"^ with a total of n vertices, find a shortest (Euclidean) 
path between a source point s and a target point t while avoiding the interior 
of the obstacles. There are two commonly studied versions of this problem. One 
is the single pair version, that asks to find a shortest path between two given 
query points. The other one is the single source version, which first constructs 
a shortest path map with respect to a source point s. After that, for any given 
query point t, a shortest path between s and t can be found based on the shortest 
path map. 

In TZ^, shortest paths are polygonal and turn only at the vertices of the 
polygonal obstacles. Sharir and Schorr 22 have developed an O(n^logn) time 
algorithm based on discrete graph searching and the visibility graph of the obsta- 
cles, where n is the number of the obstacle vertices. Various studies, e.g. [71IMT7]. 
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improved the time to quadratic in worst case. Kapoor, Maheshwari, and Mitchell 
[TU] gave an interesting 0{n + h^logn) time and 0{n) space algorithm based 
on the visibility graph approach, where h is the number of holes (obstacles) of 
the given input. This algorithm is the only algorithm known to be linear in n 
in both time and space. However, the time dependence on h is quadratic, so 
the algorithm does not perform well if h is not relatively small compared to n. 
Avoiding the visibility graph approach, Mitchell 131 developed a version of the 
continuous Dijkstra method and obtained the first subquadratic, 0(n'^/^+'^) time 
algorithm. Subsequently, based on the same technique, this result was improved 
by Hershberger and Suri [12] to O(nlogn) time and O(nlogri) space. 

In TZ^ , shortest paths among polyhedral obstacles are polygonal and turn only 
on obstacle edges or vertices. However, unlike the case in TZ^, shortest paths need 
not lie on any discrete graph. Sharir and Schorr 122 , have shown that shortest 
paths in TZ^ are geodesic, i.e. paths must enter and leave an edge at the same 
angle. Given a distinct sequence of edges, the local optimal path between two 
points can be unfolded at each edge to form a straight line, and the local optimal 
path can be uniquely identified. Nevertheless, the problem is still significantly 
harder than in TZ^. For algebraic considerations, Bajaj [5J[5] has shown that the 
algebraic complexity is exponential, since comparing the lengths of two paths 
may require exponentially many bits. Considering the combinatorial aspect of 
the problem. Canny and Reif [3] have shown that the shortest path problem in 

is NP-hard. 

An interesting special case of the shortest path problem is that in which the 
path is restricted to the surface of a single polytope. The first significant study 
of this special case in computational geometry is by Sharir and Schorr [22 . They 
gave an 0(n'^ log n) time algorithm for convex polytopes by exploiting the special 
structure of geodesic paths along the surface of a convex polytope. Mount [15] 
gave an improved algorithm for convex polytopes with running time 0{n^ logn). 
For general nonconvex polytopes, Mitchell, Mount and Papadimitriou |16] pre- 
sented an 0(n^ logn) algorithm extending the technique of Mount |T5]. See 
also [UISKHKID] for more studies for this topic. 

Our real interest is in the general case in TZ^. Several papers [6|[Tl]fT8] 
have presented polynomial time approximation algorithms for computing an 
e-approximate path, which has length at most (1 + e) times the length of the 
shortest path between two query points, where e is a real parameter defining 
the quality of the approximation. We will discuss those relevant to this paper in 
Section [TT21 

Throughout the paper we use the following notations. 

iS* denotes a set of polyhedral obstacles; 
E denotes the set of edges of S; 

Pst denotes a shortest path from s to t that avoids the interior of the obstacles; 

Qst denotes an approximate path of Pst (see next section); 

\P\ denotes the length of a given path P; 

st denotes the line segment between points s and t: 

\st\ denotes the Euclidean distance between s and t. 



1.1 Our Results 



In this paper, we address the point-to-face shortest path problem in 7?.'^: Given 
a set of polyhedral obstacles with a total of n vertices, a source point s, and a 
destination face / of some obstacle in the input, compute an e-approximate path 
from s to f that avoids the interior of the obstacles. We assume the obstacle 
surfaces are triangulated, so / is a triangle in 7?.^. 

The problem has multiple applications, including robot navigation, path 
planning and structural proteomics. For example, in structural proteomics, after 
a protein surface has been segmented and pockets identified, various descriptors 
can be associated with key points of the pocket. Given a key point s, one can 
measure the pocket depth of s by computing various distance measures from s 
to the "caps" of the pocket, the most natural of which is the shortest collision 
free distance from s to the caps (a cap is a face of the convex hull of the protein). 

We present three approximation algorithms for the point-to-face shortest 
path problem that take 0{n'^{L + \og{n / e)f / + n'^{L + \og{n / e)) / e') time, 
0{Tp-p{n) * (1/e^)) time, and 0{n'^\{n)\og{n/ e)/ + 'n'^\og{np)\og{n\ogp)) 
time, respectively, where L is the precision of the integers used, Tp^p{n) is the 
time complexity of the point-to-point shortest path algorithm used, p is the ra- 
tio of the length of the longest obstacle edge to the Euclidean distance between 
s and /, and A(n) is a very slowly-growing function related to the inverse of 
the Ackermann's function. The main contribution of this paper is in the third 
algorithm. 

1.2 Previous Work in TZ^ 

In the approximate shortest path problem, an additional real positive parameter 
e > 0, which defines the quality of the approximation, is given as part of the 
input, and the goal is to find a path between two given points s and t that 
avoids (the interior of) the obstacles and has length at most (1 + e) times the 
length of the shortest obstacle avoiding path between those two points. Such an 
approximate path is referred to as an e-approximate path or an e-approximation 
of the shortest path. In this paper, we assume e < 1. There is also a more 
general approach, where given a point s, a shortest path map with respect to s 
is constructed. In TZ^, this approach has been investigated for shortest paths on 
polytopes, as well as for shortest paths among obstacles [TT]. 

Papadimitriou [18j gave the first fully polynomial time approximation scheme 
for the general shortest path problem in TZ'^. The time complexity of the algo- 
rithm is 0(n'*(L-|-log(n/e))^/e^), where n is the complexity of the set of obstacles 
S, i.e. the number of edges, and L is the precision of the integers used, that is, 
the number of bits in the largest integer describing the coordinates of any scene 
element. 

The approach is relatively simple and thus it could be implemented in prac- 
tice. It discretizes the problem by breaking every edge into a number of small 
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Fig. 1. A shortest path and its approximation. 



segments. A visibility graph G is constructed, in which each node represents a 
segment. A Unk between two nodes is created if the two segments represented 
are visible to each other. The weight of the link is set to the Euclidean distance 
between the mid-points of the segments (see Fig. [T] for an illustration). 

To partition one edge e of S", a coordinate system is chosen such that e lies 
on the X-axis and the origin is the closest point on e to s. A sequence of points 
is added on e with coordinates {xi = eia(l + ei)*~^ : i — 1,2,3,.. .}, where a 
is the Euclidean distance from s to e and ei is a real positive parameter. As a 
result, the length of each segment is no more than ei times the distance from s 
to the segment. It is shown in [IS] that applying Dijkstra's algorithm on G gives 
a path which is at most (1 +Anei) times the length of the shortest path. Taking 
ei = e/4n gives an e-approximate path. 

Given a source point s, a real positive parameter e, and a set of polyhedral 
obstacles S in TZ^, with a total of n vertices, Har-Peled [IT] presented a technique 
to construct approximate shortest path maps in 0{Tp-p(n) * (1/e^ log 1/e)) time 
on each face and in 0{Tp-p{n)*{l/e^ log l/e)n^) time total in TZ^, where Tp-p{n) 
is the time complexity of the point-to-point shortest path algorithm used. In |11] , 
the point-to-point shortest path algorithm used is that of Clarkson [6j , which we 
will discuss later. Once the map is constructed, it takes 0(logn/e) time for each 
query, that is, given a query point t G Ti? , the length of a (1 -I- e)-approximate 
path from s to i can be reported in 0(log(n/e)) time. 

To construct a shortest path map for s on an obstacle face /, Har-Peled's 
algorithm places a set of 0(1/ log 1/e) points on /. A weighted Voronoi diagram 
is then constructed on those points. The weight of each point is the length of 
an e/%- approximate path from s to the point, obtained by any existing shortest 
path algorithm. For a point i G /, the distance function f{t) is defined as f{t) = 
Wp + \pt\, where p is the point in the weighted Voronoi diagram closest to Wp is 
the weight of p, and \pt\ is the Euclidean distance between p and t (see Fig. [3]). 
It was shown in [11] that if the points are placed "carefully", f{t) is at most 
(1 -|- e) times the length of the shortest path between s and t. 



Fig. 2. The distance function f{t) is defined as f{t) = Wp + \pt\ 



Clarkson [6] gave an 0{v?X{n) \og{n/ e) / e'^+n? \og{np) log(nlogp)) algoritlim 
for computing an e-approximation of the shortest obstacle avoiding path be- 
tween two given points s and t in 7?,'^, where p is the ratio of the length of the 
longest obstacle edge to the Euclidean distance between s and t, and A(7i) is a 
very slowly-growing function related to the inverse of the Ackermann's function. 

Let S be the set of obstacles. The general idea of Clarkson's approach is to 
construct a visibility graph V. V contains s, t, and the points on S that Pst might 
pass through. To limit the number of edges per node in y , a cone structure C is 
applied on every node x (see Fig. [S]). In a cone Ca G C, with apex x, if y is 
the closest node visible to x then an edge e between x and y is added to V. The 
weight of e is \xy\. Notice that there could be many nodes visible to x in Ca, but 
at most one edge is added. Clarkson showed that to ensure that the path found by 
running a single source shortest path algorithm on V is an e-approximate path of 
Pst, the size of the cone structure needed is 0{l/e'^^-^), where d is the dimension 
of the problem (i.e., d = 3 in T^'^). Hence, there are 0{l/e'^) edges incident to a 
node in y [5] . In the three-dimensional case, Pst could make a turn at any point 
on any edge e G E. Let a and b be the end points of an edge e E E. A point x G e 
can be expressed as a function x — f3a + {l — /3)b, where < (3 < 1. For a cone Ca 
with apex x = Pa+ {1 — (3)b, define the distance function ft,*-^°(/3) — \xy\, where y 
is the closest node in Ca visible to x. One important observation Clarkson made 
is that /i*^" is piecewise linear. Given a predefined cone structure C, each e G E 
can be divided into segments according to this piecewise linearity. The set of 
line segments formed is called the combinatorial characterization of ^ [6] , and is 
denoted as K^. According to the study of Davenport- Schinzel sequences [21], the 
size of is 0(n^A(n)/e^), where n is the number of obstacle edges, and A(n) is 
a very slowly-growing function (see above). Hence the connectivity relations of 
V can be represented by K^. 

V contains a set of carefully selected Steiner points chosen as follows. All 
endpoints of segments in are Steiner points. Additional Steiner points are 



Fig. 3. A simple cone structure with 8 cones. 



added to further divide each edge into segments no longer than -^B, where B 
is a lower bound on |Pst|- 

The algorithm for computing an e-approximate path consists of two phases. 
In the first phase, B is set to \st\, e is set to 1/2, and a simpler cone structure is 
used to compute a 1/2- approximate path. This provides a better lower bound 
on \Pst \ for the second phase, which gives the final approximation. 

2 Approximate Point-to-Face Shortest Paths 

In this section we present three algorithms for finding approximate point-to-face 
shortest paths among polyhedral obstacles in TV'. 

2.1 A General Approach 

The main idea used is to place (a grid of) Steiner points on the target face / such 
that there is at least one point close enough to t to give us a good approximation. 

Let B and D denote a lower bound and an upper bound on the shortest path 
between s and /, respectively. 

Lemma[I]below is an application of Lemma 2.10 of [llj . which proved a more 
general claim. 

Lemma 1. Let x be any point on the target face f . Let h be the point on f closest 
to the source s with respect to the Euclidean distance. We have \Psh\ £ 2[Ps2.|. 

Proof. Refer to Fig. SI Let s' be the projection of s on the plane supporting /. 
If s' £ /, then h = s' , otherwise h lies on the boundary of /. It is obvious that 
\hx\ < \sx\ < \Psx\. We have \Psh\ < \Psx\ + \hx\ < 2\P,^\. □ 



face 



Fig. 4. \Psh\ is no more than 2|P5a:|. 



Let t G / be the end point of the shortest path between s and /. Hence the 
shortest path between s and / can also be expressed as Pgt ■ By Lemnia[Tl we have 
\Psh\ < 2|Psf|. By computing an e-approximate path between s and h, i.e. Qsh-, 
we can obtain the lower bound and upper bound on \Pst\ as B = |(5s;i|/(2(l + e)) 
and D \Qsh\- 

2.2 The First Algorithm 

Our first algorithm for finding approximate point-to-face shortest paths extends 
Papadimitriou's algorithm for the point-to-point version. We proceed as follows. 
The lower bound B and the upper bound D of the shortest path can be computed 
by identifying the point h and computing Qsh using Papadimitriou's algorithm. 
Obviously \ht\ < Pgt < D, which implies t must be within distance D oi h (see 



We define a sample grid as a uniform grid of unit length eB, which is applied 
on the target face / within distance D of h. The total number of sample points 



Lemma 2. There exists one sample point t' that gives a Se-approximate path. 
Proof. Let t' be the sample point closest to t. We have \tt'\ < eB < e\Pst\ and 



Fig. ED. 



is 0{l/e^). 



\Pst' 

\Qst' 



< \Pst\ + \tt'\ 

< {l + e)\Pst'\ 



< (l + e)(|P,i| + |tt'|) 

< {1 + e){\Pst\ + e\Pst\) 

< (l+3e)|P,t| 



where Qst' is the approximating path and Pgt is the shortest path. 



□ 



Fig. 5. The point t must be witliin D distance of h. 



In Papadimitriou's algorithm, each edge of S is split into at most 0{{L + 
log(l/ei))/ei) segments [TH]- The number of visibility graph edges resulting from 
this subdivision of edges of S is 0((n((L + log(l/ei))/ei))^). With ei = e/4n we 
have 0{n'^{L + \og{n/e)f/e^) edges. We also have 0{n^{L + log{n/e))/e)*{l/e^)) 
graph edges corresponding to visibility edges between the segments on edges of 
S and the sample points on /. The number of vertices of the graph is 0{n?{L + 
log(n/e))/e + 1/e^). Dijkstra's algorithm takes 0{N log N + M) time, where M 
is the number of graph edges and N is the number of graph nodes. With M = 
0{n'^{L+\og{n/e)f/e'^W{L+\og{n/e))/e^) and N = n'^{L+\og{n/ e))/ e+l/ e"^), 
Dijkstra's algorithm takes 0{n^{L + log(n/e))^/e^ + v?{L + log(n/e))/e^) time 
to compute a "it-approximate path between s and / by Lemma [21 

In order to have an e- approximate path between s and /, we simply use a 
new parameter e" = e/3 to construct the visibility graph. The new parameter 
does not affect the complexity of the algorithm. 

Theorem 1. Given a point s, a face f, a set S of obstacles, with a total of 
n vertices, and a real positive parameter e, an e-approximation of the shortest 
obstacle-avoiding path from s to f can be computed in 0(ri^(L + log(n/e))^/e^ + 
n^(L + log(n/e))/e'^) time, where n is the complexity of S and L is the number 
of bits in the largest integer describing the coordinates of any scene element. 

2.3 The Second Algorithm 

Our second algorithm builds upon Har-Peled's point-to-point approximate short- 
est path algorithm. Har-Peled' algorithm constructs an approximate shortest 
path map on a given face / with respect to a source s in 0{Tp-p{n)*{l / e'^ log 1/e)) 
time, where Tp-p{n) is the complexity of the point-to-point shortest path algo- 
rithm used. As mentioned previously, a full map on / is unnecessary as t must 
locate within D radius of h. Hence only a partial map within D radius of h is 



sufficient to capture a good approximation. There is no actual query to be per- 
formed, since the sample points placed on / to compute the weighted Voronoi 
diagram already serve the purpose as a sample grid. The number of sample 
points placed within D radius of h is 0(l/e^) ^11,. It takes 0{Tp-p{n) * (1/e^)) 
time to construct the partial map. If Clarkson's algorithm is used for Tp_p, we 
have: 

Theorem 2. Given a query point s, a face f and a real positive parameter 
e, an e-approximation of the shortest obstacle-avoiding path from s to f can 
be computed in 0((n^A(n) log(n/e)/e^ + log(np) log(n logp)) * (1/e^)) time , 
where p is the ratio of the length of the longest obstacle edge to the Euclidean 
distance between s and t, and X{n) is a very slowly- growing function related to 
the inverse of the Ackermann's function. 

This second algorithm is very similar to the first one. Both algorithms need 
a sample grid of size 0(l/e^) and execute a point-to-point shortest path algo- 
rithm on each point, which means we have to run a point-to-point shortest path 
algorithm 0(l/e^) times. In the next subsection, we show how to obtain a point- 
to-face shortest path algorithm which has the same asymptotic complexity as 
Clarkson's point-to-point version. 

2.4 The Third Algorithm 

Theorem 3. Given a query point s, a face f , a set of obstacles S and a real 
positive parameter e, we can find an e-approximate path of the shortest obstacle- 
avoiding path from s to f in time 0{n^X{n) \og{n/e)/e^ + n'^ log(np) log(n log p)), 
where p is the ratio of the length of the longest obstacle edge to the Euclidean 
distance between s and t, and X{n) is a very slowly- growing function related to 
the inverse of the Ackermann's function. 

Notice that in Theorem[3]the point-to-face approximation has the same complex- 
ity as Clarkson's point-to-point approximation 6. . Indeed, our solution builds 
upon Clarkson's point-to-point solution. 

We will use a different approach to prove Theorem [3l Instead of building a 
sample grid on / first, we will construct a visibility graph first and then decide the 
additional Steiner points needed on / based on the visibility graph constructed. 
All additional Steiner points on / can be added directly to the visibility graph. 
Hence, we only need to execute Dijkstra's algorithms once for all Steiner points 
placed on /. We will explain how additional Steiner points and the corresponding 
visibility graph edges are produced, and give a count on the total number of 
Steiner points. 

Clarkson's algorithm consists of two phases. In the first phase, a coarse ap- 
proximation of the shortest path between the source and destination is computed, 
which gives a better lower bound for the second phase. We set h, which is the 
closest point on / to s in Euclidean distance, as the destination point and find a 
5 — approximation of the shortest path between s and h using the first phase of 



Fig. 6. Pst cannot travel beyond the sphere centered at h with radius 2D. 



Clarkson's approach. Let t G / be the end point of the shortest path between s 
and /. Set the length of the ^ — approximate path as the upper bound D of Pst 
and let B — -0/2(1 + 1/2) as the lower bound, following Lemma. [TJ Obviously, 
Pst cannot travel beyond the sphere centered at h with radius 2D (see Fig. [S]) . In 
the second phase, following Clarkson's approach, we only need to partition ob- 
stacle edges within the sphere into segments of length eB/An and then apply the 
cone structure on each Steiner point created. With this, we have the following 
lemma. 

Lemma 3. // 1 is on an edge e G / then there exists a Steiner point t' £ e such 
that IQst'l < (l + 5/4e)|Pst|. 

Proof. On e, choose the Steiner point closest to t as t' . See Fig. [7] for an illustra- 
tion. Since each segment is no longer than ^|Psf|, \tt'\ < ^jPstj. By triangle 
inequality, we have |Pst'| < |Pst| -I- \tt'\ and 

\Qst' \ < (1 + e)|Pst'| {from Clarkson s Steiner point placement) 

<(l+6)(|PstK|tt'|) 
<(l+6)|Pst| + (l + e)(^|Pst|) 

< (1 -I- 5/4e)|Pst| [since e < 1 and n>l) 

□ 
□ 

To handle the case when t is in the interior of / we need to add two sets of 
Steiner points. We do this as follows. For the first set, for each obstacle edge e, 



Fig. 7. The Steiner point t' gives a good approximation for t. 



we add point u G e as a new Steiner point, if u satisfies: 

1) its projection u' in the plane containing / is in the interior of /, and 

2) u is visible from u', and 

3) uu' is tangent to some obstacle, which does not contain e. 

We apply the cone structure (following Clarkson's algorithm) on u and add the 
corresponding edges to the visibility graph. Since there are at most such 
points and the number of Steiner points in the original visibility graph is J7(n^), 
the complexity of the visibility graph is unchanged. 

To obtain the second set of Steiner points, we proceed as follows. For each 
existing Steiner point v on an obstacle edge, let v' be the projection of v in the 
plane containing /. We add v' and edge v'v to the visibility graph if v' is in 
the interior of / and v is visible from v' . Note that we do not apply the cone 
structure on v' . Observe that, we at most double the number of Steiner points 
and each additional Steiner point introduces exactly one edge. The complexity 
of the visibility graph remains the same. Wc will discuss later how to efficiently 
determine the visibility between a point and its projection on a plane. 

Lemma 4. Let the contact point of Pgt on the last obstacle edge before reaching 
f be ci, that is, cit is the last segment of Pst- If t is an interior point of f then 
c\t must be perpendicular to f . 

Proof. We make the proof by contradiction. Suppose t is in the interior of / and 
Cit is not perpendicular to /. Wc will construct a new path from ci to some 
point t' G / such that the new path is shorter than \cit\. 

Let c'^ be the projection of ci on the plane tt containing /. Consider the plane 
formed by ci, c'^ and t. Let t' = t. We shift t' along the segment tc'i until either 
c[ is reached or the segment cit' intersects some edge e at a point C2. We have 
\cit'\ < \cit\, a contradiction. 

□ 



Cl' t' t 

Fig. 8. Illustration of (the proof of) Lemma H) 

Lemma 5. If t is an interior point of f , in the modified visibility graph there 
exists at least one Steiner point t' on f and such that \Qst'\ ^ (1 + 5/4e)|Psj|. 

Proof. By Lemma^ the last segment cit of Pst is perpendicular to /, since t is an 
interior point of /. Assume ci is not a Steiner point, since otherwise we are done. 
On the edge containing ci, let di be one of the two Steiner points neighboring 
Cl, specifically, one visible from its projection on /. Such a Steiner point exists 
from our placement of Steiner points (set one). We have \cidi\ < ^|Ps(|. Let the 
projection of c?i on f he t' . Use an alternative path from s to t' to approximate 
Pst (See Fig. |9] for an illustration). We have 

\Qst' \ = \Qsd,\ + \dit'\ 

< \PsdA+^\Pst\ + \dit'\ 

< \PsdA+^\Pst\ + \cit\ + Icidil 

< \Pst\ + e\Pst\ + -^\Pst\ 

An 

< (l + 5/4e)|Pst| 

□ 

Following Lemma [3] and Lemma EJ the new visibility graph gives us a 5/ 4e- 
approximate path from s to f. To conclude the proof of Theorem [3] we need to 
show how to find the additional Steiner points and visibility edges within the 
given time bound. 

Visibility Computation We now show how to compute the additional Steiner 
points on obstacle edges as well as on the target obstacle face /. Since addition of 
Steiner points does not change the complexity of the visibility graph, the cost of 
the shortest path computation is 0{n^X{n) log(n/e)/e'* + n^ log(np) log(nlogp)). 
We would like to keep the computation of additional Steiner points and edges 
that capture the visibility between Steiner points on obstacle edges and their 
projections on / within this time bound. 
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Fig. 9. Approximation of P^t by an alternative path. 



A naive approach to compute the visibility between two points u and u' 
would require 0{n) time by checking whether the line segment uu' intersects 
any obstacle. Over all points, this will exceed the time bound above. 

We can try to reduce the time to compute the visible pairs by using a ray 
shooting data structure. As mentioned in [53], the general ray shooting problem 
in three dimensions is still far from being fully solved. However, our problem is 
a special case of the ray shooting problem. We need to determine the visibility 
between a point and its orthogonal projection on the target plane tt supporting 
/, instead of the visibility between any two arbitrary points, i.e. the direction for 
ray shooting is always perpendicular to n. On each obstacle edge e, all Steiner 
points and their projections on w are coplanar. Let tt' be the plane passing 
through e and perpendicular to tt. By computing G = tt' P| 5* for each edge of 
S (see Fig. [T0|) . we can transform the three dimensional visibility problem into 
0{n) subproblems in dimension two, where n is the number of obstacle edges. 

Let edge e' be the projection of an edge e onto tt. Notice that all line segments 
between Steiner points on e and their projections on tt (that is, on e') are parallel. 
We can perform a plane sweeping on G with a horizontal line and maintain a 
data structure which stores edges intersecting with the sweep line at a certain 
time instance. The data structure is updated when the sweep line passes a vertex 
and each insertion and deletion precess requires O(logn) time [19j . During this 
process, we could partition e into 0{n) open segments, such that all points on 
each segment have the same visibility answer with respect to their corresponding 
projections on e'. Furthermore, the end points of segments that are visible from 
their projection on / are the new Steiner points to be added on obstacle edges. 
The total time of adding Steiner points is 0{n^ \ogn+N), where N is the number 
of Steiner points in the original visibility graph. 

We can simplify this process as follows. Since the order of the edges inter- 
secting with the sweep line is not important, we can avoid the plane sweep 
procedure altogether while still partitioning each edge into 0{n) segments in 



Fig. 10. The sweep line intercepts edges e, 1, 2, 3, 4, and q 



O(nlogn) time. Recall that we assume the obstacles are triangulated. For each 
triangle A, find the intersection segment r = tt' P| Z\. If at least one vertex of 
r is inside of the quadrilateral formed by the endpoints of e and e', we project 
r on e. The projection of r on e is the region in which the visibility of Steiner 
points to their counterparts are blocked by r. Hence we can divide e into 0(n) 
segments in 0(n log n) time. Again, the total time is 0(n^ logn + N). 

Note that the cost of visibility computation is dominated by the cost of the 
shortest path computation, i.e. the overall time to approximate the shortest path 
from s to / is 0{v?\{n) log(n/e)/e^ + log(np) \og{n\ogp)). 

This concludes the proof of Theorem [3l 

3 Conclusions 

In this paper we discussed three point-to-face approximate shortest path al- 
gorithms in Ti? . It is interesting to notice that the point-to- face shortest path 
must end on / within a certain range from h, where h is the point on / clos- 
est to s in Euclidean distance. That range can be obtained by applying known 
point-to-point shortest path algorithms between s and h. After placing a sample 
grid near /i, we proved we can find an e-approximate path between s and / 
in 0{n'^{L + log(n/e))^/e^ + + log(n/e))/e^) time by extending Papadim- 
itriou's algorithm or in 0((n^A(n) log(n/e)/e'* -I- n? \og{np) log(nlog p)) * (1/e^)) 
time by extending Har-Peled's algorithm. However, this approach still requires 
to execute a point-to-point shortest path algorithm 0(l/e^) times. Finally we 
showed that Clarkson's point-to-point shortest path approach can be extended 
to solve the problem by adding additional Steiner points directly to the visibility 
graph, without changing the asymptotic complexity of the algorithm, resulting 
in an 0(n^A(n) log(n/e)/e^ -I- in? \og{np) log(nlogp)) time algorithm for finding 
an e-approximate path between a source point s and an obstacle face /. 
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